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OHMIC  EFFECTS  IN  QUASIOPTICAL  RESONATORS 


I.  Introduction 

The  quasioptical  gyrotron  (QOG)  has  been  under  investigation  at  NRL  for  several  years 
as  a  potential  source  of  high-power  radiation  at  frequencies  above  100  GHz  for  electron 
cyclotron  resonance  heating  (ECRH)  of  fusion  plasmas.1’ 2'  3’  4  An  integral  part  of  the  QOG 
is  the  resonator,  which  traps  some  of  the  radiation  emitted  by  the  electron  beam,  allowing 
the  radio  frequency  (rf )  fields  to  grow  to  the  large  values  necessary  for  efficient  extraction  of 
energy  from  the  beam.  I  lie  QOG  resonator  surfers  losses  due  to  diffraction  of  the  radiation 
around  the  outer  edge  of  the  resonator  mirrors  as  well  as  the  finite  conductivity  of  the 
mirror  surfaces.  Clearly,  it  is  imperative  to  understand  each  of  these  loss  mechanisms  in 
order  to  understand  the  operation  of  the  QOG  experiments.  This  note  addresses  the  ohmic 
losses  in  the  resonator,  deriving  expressions  for  the  electric  field  in  the  resonator,  the  ohmic 
quality  factor  (Qn)<  the  ohmic  heating  density  (pa).  the  total  ohmic  power  (Pq)  dissipated 
in  the  resonator  mirrors,  and  the  temperature  rise  of  the  mirrors.  The  derived  formulae  are 
then  applied  to  the  current  design  of  the  cavity  to  be  utilized  in  the  NRL  induced  resonance 
electron  cyclotron  (IREC)  maser  experiment. 


II.  Cavity  Fields 

The  resonator  considered  here  consists  of  the  Fabry-Perot-type  open  resonator  shown 
in  Fig.  1.  The  spherical  cavity  mirrors  form  an  azimuthally  symmetric  resonator  about 
the  cavity  ( y' )  axis,  which  is  offset  by  an  angle  6  from  the  y-axis  as  shown  in  Fig.  2.  The 
angle  0  is  used  in  the  QOG  literature,5  however,  in  the  IREC  maser  literature6  the  angle 
between  the  resonator  (y')  and  the  electron  beam  (z)  axes  is  typically  used.  The  mode 
structure  and  stability  of  this  type  of  resonator  is  discussed  by  Yariv7  who  finds  that  the 
electric  field  of  the  transverse  electric  and  magnetic  (TEMm>n  /)  modes  of  the  travelling 
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DEPRESSED 


Figure  1:  Schematic  diagram  of  the  quasioptical  gyrotron  experiment. 


GiienUtion  of  the  resonator  (t/')-axis  with  respect  to  the  electron  beam  (z)- axis 


wave  in  the  cavity  may  be  expressed  as 


E  =  xEt 


where  the  physical  electric  field  is  obtained  by  taking  the  real  part  of  E,  and 


r»t  ,  ,  ,  _  E^J  WfJ  (  S/IX  \  „  (  -fU  \  /_  (fl+fll  \ 

)  2  ”U-to')jeXPi  w?(sr')  / 


Eq,i  Wq,i  jr  (  y2x  \  /  V2z'  \ 

2  MW)  \MV))  \m(y,))e 

xexpl-.  L 


(m  +  n  +  1)  arctan 


(911 


where  Hm  is  a  Hermite  polynomial  of  order  m,  wQ  is  the  radiation  beam  waist, 

( y )  =  <,  (i  +  y2/zl) , 

#™(y)  =  y(i  +  ^/y2), 

7ru'o,/  d  / 1  +  g\  ! 


\  2  vi -W  ’ 

72^  is  the  radius  of  curvature  of  the  radiation  wavefront,  fc  is  the  wave  number,  A  is  the 
wavelength  and  u;  is  the  angular  frequency  of  the  radiation,  d  is  the  resonator  mirror 
separation  and  g  —  1  —  d/7?c,  where  Rc  is  the  mirror  radius  of  curvature,  which  is  assumed 
equal  to  the  curvature  of  the  radiation  wavefront  at  the  mirror.  The  electric  field  of  the 
standing  wave  in  the  resonator  is  obtained  by  adding  the  two  oppositely  directed  travelling 
waves,  yielding 


My')/  I  W(y') 


(x2  +  z'7) 


x  exp  j uit  +  (m  +  n  +  1)  arctan  ^ — 

(  (x2  +  2 '*)  /tt] 

x  cos  ki  y  +  — - - — -r - , 

V  2  Rw(y')  2 


where  the  fa/2  term  in  the  cosine  is  necessary  to  satisfy  the  boundary  conditions  at  the 
resonator  mirrors. 
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Of  particular  interest  is  the  electric  field  at  the  electron  beam,  which  is  located  very  near 
the  RF  beam  waist  (i.e.  y'  z0 ).  In  this  case  «  0,  tuj(y')  ss  w0,  and  arctan(y'/2<>)  %  0 
yielding 


F3 

m,n,l 


(■ x,y',z',t )  =  EajHm 


x  exp  {luJit}  cos 


(7) 


For  the  TEMo.o,/  mode  this  becomes 


=  E0,i  exp 


|  exp  cos 


(8) 


The  peak  electric  field  of  the  TEMo,o,i  mode  along  the  resonator  (y1-)  axis  is  (in  MKS  units) 


vpeak 


Jo,l 


Wo,l 

ivi(y') 


_  m0cu 7  w0j  {  sin2(i/>)  \ 

e/?r4  wi(y')  \(1  -0 ||Cos(^))2J  V  sa 

where  F  is  the  peak  normalized  wave  amplitude  commonly  used  in  the  gyrotron  literature,8,  9 
m0  and  e  are  the  rest  mass  and  charge  (magnitude)  of  an  electron,  7  is  the  relativistic  fac¬ 
tor,  and  Px  and  are  the  electrons  velocity  perpendicular  and  parallel  to  the  magnetic 
field  normalized  to  c,  the  speed  of  light.  The  parameter  s  is  the  harmonic  number  of  the 
beam-wave  interaction  for  the  case  of  the  QOG  and  should  be  set  to  unity  for  the  IREC 
maser  resonator.  The  angle  xp  is  the  angle  between  the  resonator  and  electron  beam  axes, 
and  should  be  set  to  7r/2  for  the  QOG. 


III.  Ohmic  Power  Losses 

To  calculate  the  ohmic  power  lost  in  the  resonator  mirrors,  consider  a  plane  wave 
incident  on  a  semi-infinite  planar  conductor  as  shown  in  Fig.  3.  Assuming  that  all  field 
quantities  S  vary  in  time  as 

5(x,  t)  =  S(x)e,u"  (10) 
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Figure  3:  A  plane  wave  incident  on  a  semi-infinite  planar  conductor. 


and 


B 

D 

P 

J 

Maxwell’s  equations  become 

V  x  E 

V  x  H 

V  D 

V  B 

Using  a  vector  identity  we  obtain 


fill  (11) 

eE  (12) 

0  (13) 

<rE,  (14) 

-iw  B  (15) 

(iu;e-|-a)E  (16) 

0  (17) 

0.  (18) 


=  V  (V  •  E)  -  V2E 


V  x  V  x  E 
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(19) 


=  — V2E 


=  — iw^tV  x  H, 


from  which  the  wave  equation  follows: 


V2E  =  iufj.  ( iuit  +  <r)  E. 


Similarly,  one  may  derive  the  wave  equation  for  H 


V2H  =  iufj.  (iue  +  a)  H. 

Region  1  (see  Fig.  3)  is  free  space,  so  that  a  =  0,  leading  to 


V2Ei  =  —u2  HotoEi. 


Picking  the  orientation  of  the  electric  field 


Ei  =  xEi{x), 


the  solution  to  Eq.(24)  is 


where 


Ei  =x{E+e~iky  +  E;eiky} 


Ic  —  LO  HqCq 


In  region  2  we  assume  ue  -C  cr  (i.e.  the  displacement  current  is  negligible)  yielding  the 


wave  equation 


where 


V2E2  =  iujfiaE/2  =  t2  E2 


where  6  is  the  skin  depth  of  the  conductor.  The  solution  to  Eq.(28)  is  then 

E2  =  x{Ele~Ty  +  E;eTy) 

=  xEf  e-TV 
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(30) 

(31) 


(41) 


The  time  averaged  power  per  unit  area  lost  in  the  conductor  (p n)  is  given  by  the  real 
part  of  the  Poynting  vector  (N)  evaluated  at  x  =  0; 

ft>  =  »f(N)|„0}  (42) 


where 


<N)Uo  = 


2^2  x  H; 


=  y~-{l  +  i)\E  2+|2 


=  y 


4(1  +i)E+ 
rr6  [(Z0  +  2 /<t6)2  +  Z 2 


For  reasonably  good  conductors  and  frequencies  not  tc  >  high,  the  factor 

2  /2u?/i0\ 1/2 


is  small  compared  to  Z0  (=  377f l)  and  the  heating  density  becomes 

(2EtV 

*  *  W  '  <47) 

In  addition  to  being  applicable  to  the  case  of  a  plane  conductor,  this  equation  is  a  good 
estimate  of  the  heating  density  for  a  curved  conductor  when  the  radius  of  curvature  is 
much  larger  than  the  wavelength  of  the  radiation.  Using  Eq.(2)  for  a  TEM0,o,(  mode  at  the 


mirror  we  have 


2E t  =  «P  {-rV«>?  (rf/2)}  • 


1  '  w,(dl 2)  r  l  '  '  1  ' 

Therefore,  the  power  lost  per  mirror  due  to  ohnne  heating  in  a  mirror  of  radius  a  at  position 


=  I  Pn^r 
Jo 


*Ehwl,i 


4(7  bZl 


1  —  exp 


wf  (d/2) 
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IV.  Stored  Energy 


In  vacuum,  the  time-averaged  energy  density  is 

u  =  j  (E  •  D*  -f  B  •  H')  =  •  D’  (50) 

The  total  stored  energy  in  the  resonator  is  twice  the  stored  energy  due  to  the  wave  travelling 
in  one  direction  given  by  Eq.(2).  For  a  TEM0,o,/  mode  this  yields 

tjwljEl,  rH 2  [°  2*r  dr 


^stored  — 


toKjKi  r'1  d  \a  { 
4  J-d.li  ^  Jo  i 

*£o Ehwljd  [  ( 

=  - § -  i1  ”  eXP  ( 


wf  (y) 

'  -2a2 

wl(d/ 2) 


V.  Ohmic  Q 

The  ohmic  Q  of  a  cavity  is  defined  as 

Qn  = 


OJ  Estored 


Noting  that  power  is  lost  in  two  mirrors,  and  using  Eqs.(49)  and  (51)  we  obtain 

Qn  =  \  {*fVoe)1/2 

where  /  is  the  frequency  of  the  radiation. 


VI.  Resonator  Heating 

In  this  section  we  wish  to  calculate  the  time  dependent  rise  in  temperature  of  the 
resonator  mirrors.  To  do  so,  we  solve  the  one-dimensional  heat  equation  with  appropriate 
boundary  conditions11: 

A 

ut{x,t)  -  auri(x,t)  =  —  (54) 

pc? 

u(x,.  0)  =  u0  (55) 

ur(0, t)  =  0  (56) 

U:{LJ)  =  0  (57) 

A  =  Tb{x)  (58) 
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where  u(x,t)  is  the  temperature  of  the  mirror  at  position  x  and  time  t,  the  subscripts  of  u 
denote  partial  derivatives  of  u,  the  mirror  is  located  between  x  =  0  and  x  =  Z,,  T  is  the 
flux  incident  on  the  mirror,  and  a  =  k/  (pc,)  where  p,  Cp,  and  k  are  the  mass  density,  heat 
capacity,  and  thermal  conductivity  of  the  mirror  material. 

The  first  step  in  solving  this  problem  is  to  set  T  =  A  =  0  and  separate  variables 
to  find  the  eigenvalues  and  eigenvectors  of  the  related  homogeneous  problem.  Setting 
u  =  X(x)T(t)  and  denoting  derivatives  by  primes  gives 


Therefore 


X"  +  \2X  =  0 

(60) 

X\0)  =  X\L)  =  0 

(61) 

which  has  solutions 


Xn{x)  =  cos(Anx)  An  =  n  =  0,1,2,.... 
Now  let  A  /  0  and  try  the  solution 

u{x,t)  =  ^2Tn(t)Xn{x)  . 

n 

Multiplying  by  u(x,t)  and  integrating  yields 
To (0  =  i  f0Lu(x,t)dx 

Tm(t)  ~  \  f0L  u(x,t)  cos  (mf^)  dx  m=  1,2,... 

Next  we  expand 

A/(pCj,)  =  £qn(t)Xn(x)  . 

n 

Inserting  Eq.(58)  gives 

\qo=^l 
\  qm  =  m=  i’2”-- 

Therefore  the  differential  equation  reduces  to 


(rn  +  aA2„Tn  -  qn)  =  0 


(62) 

(63) 

(64) 

(65) 

(66) 

(07) 
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Table  1:  Parameters  for  the  quasioptical  IREC  maser  experiment. 


Mode 

TEM0,o,i 

Output  power 

15  MW 

Frequency 

280  GHz 

Wavelength 

1.1  mm 

Magnetic  field 

6  T 

Electron  beam  energy 

500  keV 

Electron  beam  current 

200  A 

F 

0.2 

a  =  /?±//5|  | 

0.5 

o 

O 

C4 

12 

which  leads  to 


u0  +  got 


n  =  0 


[eXP(_AnQ*)  -  l\  U=  1,2,... 


(68) 


\za 

x  n 

Thus  the  solution  for  the  mirror  temperature  is 

,  .  Ft  2TL  ^  l  / tittx\ 

u(I’  cos  (t  J 


1  -  exp 


— n27r2a< 

Z2 


(69) 


VII.  The  IREC  Maser  Resonator 

In  this  section  the  formulae  derived  will  be  applied  to  the  current  design  of  the  resonator 
for  the  quasioptical  IREC  maser  experiment.  Representative  experimental  parameters  are 
given  in  Table  1.  The  theory  necessary  to  optimize  an  IREC  maser  design  is  complicated 
and  will  be  addressed  in  future  publications,  however,  let  us  pick  F  =  0.2  for  this  example. 
This  value  is  approximately  twice  the  F  value  for  optimum  efficiency  operation  in  a  QOG 
and  is  appropriate  for  the  IREC  maser  since  here  the  beam  electrons  interact  only  with 
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the  RF  fields  propagating  in  the  direction  of  the  electron  beam  travel.  Also  picked  are 
values  for  the  angle  between  the  resonator  and  electron  axes  ip  =  20°  and  the  normal¬ 
ized  interaction  length  p  =  12.  Specifying  the  peak  heating  density  at  the  mirror  to  be 
150  kW /cm2  determines  the  resonator  mirror  separation  (81  cm)  and  radius  of  curvature 
(52  cm).  The  output  coupling  of  14%  is  picked  to  balance  the  input  and  output  power 
levels  and  determines  the  mirror  radius  of  3.0  cm.  Assuming  values  for  the  conductivity 
(3.6  x  107  siemens/m)  and  the  skin  depth  (1.6  x  10-7  m),  appropriate  for  copper,  Eq.(49) 
may  be  used  to  calculate  the  power  lost  per  mirror  (780  kW). 

With  the  resonator  parameters  specified,  it  is  a  simple  matter  to  apply  the  derived 
formulae  to  obtain  the  peak  electric  field  at  the  center  of  the  resonator  (Eq.(9))  E0  = 
1.1  MV/cm.  Similarly,  the  energy  stored  in  the  resonator  (Eq.(51))  is  Ea tored  =  2.4  J,  and 
the  ohmic  quality  factor  (Eq.(53))  is  Qq  =  2.5  x  106.  The  resonator  heating  is  calculated 
assuming  that  the  incident  flux  is  T  —  150  kW/cm2,  the  peak  heating  density.  Eq.(69)  then 
gives  a  temperature  rise  at  the  surface  of  the  mirror  of  approximately  150°C  for  a  2  /isec 
pulse.  Standard  values  for  copper  of  p  =  8.93  g/cm3,  cp  =  0.385  J/(g°K),  k  =  3.5  W/cm°K, 
a  =  1.02  cm2/sec,  and  a  mirror  thickness  of  L  =  1  cm  were  used  in  this  calculation.  It 
should  be  noted  that  the  number  of  terms  necessary  for  the  convergence  of  the  sum  in 
Eq.(69)  is  approximately  1000.  Figure  4  shows  that  at  the  end  of  2  psec  the  temperature 
rise  in  the  mirror  is  confined  to  very  near  its  surface.  The  parameters  calculated  for  the 
IREC  maser  resonator  are  summarized  in  Table  2. 
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Table  2:  Calculated  parameters  for  the  quasioptical  IREC  maser  experiment. 


Mirror  separation  ( d ) 

81  cm 

Mirror  radius  of  curvature 

52  cm 

Output  coupling 

14  % 

Mirror  radius 

3.0  cm 

E0 

1.1  MV/cm 

W0 

8.5  mm 

w(d/ 2) 

1.8  cm 

21  cm 

Quality  factor  (diffraction,  Qa iff) 

66,000 

Quality  factor  (ohmic,  Qn) 

2.5  x  106 

Pn  P«k 

150  kW/cm 

Ohmic  power  lost  (per  mirror) 

780  kW 

13 


150 


Figure  4:  Temperature  rise  in  the  resonator  mirror  calculated  from  a  one-dimensional  model 
for  an  incident  power  flux  ( } r)  of  150  kW/cm2  and  a  2  /xsec  pulse  length. 
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